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FORWARD 


This document was prepared by personnel at Continuum, Inc. for NASA-MSFC under 
Contract NAS8-35506. This document is an interim report on the first year of work 
under thiB contract. 


1. INTRODUCTION 


1.1 B ackg r ound 

The highly non-uniform flow around the LOX posts in the SSME powerhead has 
contributed to a long history of failures of the posts. Both pressure and heating loads 
have caused problems which have resulted in undesirable, but necessary, design 
modifications such as the use of LOX post flow shields. The geometric complexity of the 
LOX post flowfield is enormous; 600 posts are fed by the five hot gas discharge ducts 
from the HPFTP and HPOTP. The posts are fluted to modify the structure of the trailing 
vortices and are shielded by plates covering pairs of posts in the outermost row. Hot 
gase6 flow along the sides of the injector elements into entry ports which conduct the 
flow through an annulus into the main, combustion chamber. The region of posts which 
are subjected to extreme environments is contained within the region bounded by the 
exits of the hot gas transfer ducts, the bottom of the oxidizer manifold, and the space 
above the secondary plate. The hydrogen cavity flow between the primary and secondary 
plates does not cause severe environments arid, is not considered further. 

Assuming that the flow from the HGM is symmetric about a plane through the center 
transfer tube, one-half of the region could be modeled at one time. Even the half-plane 
flow would be too complex to provide a direct numerical solution to the flow field of 
interest. Continuum has been contracted to address this problem by a phased effort 
which first models the flow around a single and small clusters (2-10) of posts, second 
models the velocity field in the cross-flow plane, and third models the entire flow region 
with a 3-dimensional network-type model. However, the contract has been modified to 
include a full 3-D numerical solution of the flow field in the high pressure fuel turbopump 
turnaround duct (TAD), hot gas manifold (HGM) and transfer tubes. The results of this 
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effort will be used to define boundary conditions for LOX post analyses* 

1.2 Objectives 

The following sections discuss the work performed under Phase I of the contract as 
modified to reflect the TAD/HGM/Duct analysis. These sections include a presentation 
of Continuum's laminar and turbulent boundary layer development in support of the LOX 
post study and the results of the 3-D HGM analysis# 


2. SUMMARY 

Continuum has developed shear stress wall functions which will permit viscous analyses 
without requiring excessive numbers of computational grid points* These wall functions, 
laminar and turbulent, have been compared to standard BlasiuB solutions and are directly 
applicable to the cylinder-in-crossflow class of problems to which the LOX post problem 
belongs* The results of this work are presented herein. 

Continuum has also performed a full 3-D fluid flow analysis of the HPFTP exhaust 
system which consists of the turnaround duct, 2^duct hot gas manifold and the " Version 
B" transfer ducts. The results of this analysis are presented in this report. 


3. HIGH REYNOLDS NUMBER CROSSFLOW ABOUT A CYLINDER 
3,1 Introduction 

In order to accurately account for pressure loading and heating to cylinders in crossflow, 
like for the LOX posts in the powerhead, a detailed flowfield prediction and suitable wall 
boundary conditions are required. It is impractical to resolve the flowfield in the vicinity 
of the wall with enough grid points to accurately calculate either the wall friction or 
heating, therefore a special wall treatment is required. "Wall functions" are commonly 
used to provide the required boundary conditions; however, care must be exercised in 
order not to make the wall functions too empirical. The end result is to predict 
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frictional losses from the detailed flow vectors, not from the mean channel flow. The 
following procedure was developod and tested with geometrically simple problems to 
provide the necessary CFD tools for powerhead analysis. 


3.2 Laminar Boundary Layers 

Laminar boundary layers on a flat plate were analyzed with Continuum's VAST code for 
constant density and temperature (hence pressure uIbo) conditions. Figure 3-1 shows the 
results of these analyses compared to the Blasius solution. A mean lateral velocity was 
used as a boundary condition on the free-stream side of the computational region. The 
excellent agreement of the VAST solution with the Blattius solution suggests that no 
significant artificial viscosity effects are present in the solution for this case. Notice 
that only 15 lateral node points were used for this calculation. Identical results were 
obtained for a case run with 11 nodes, only one of which was initially in the boundary 
layer. When 7 nodes were used with a step velocity profile a solution was generated, but 
the accuracy of this solution w.,s reduced. These results are considered acceptable, and 
the use of at least one node in the initial boundary layer is reasonable. 


3.3 Turbulent Boundary Layer 

For turbulent wall flows, including fully developed pipe flows, for a smooth wall the 
following empirical velocity profiles are valid. 


~ 4 _ 
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u + = 5.5 + 2 . 5 ln{y + ) for y + > 30 


(3-1) 


u + * -3.05 + 5.0 ln(y + ) for 5 < y + < 30 


(3-2) 


+ + . + / c 

u = y for y < 5 


(3-3) 


where 

+ , * 
u - u/U 

U* - <T o /p) 0,S 

/ , , 0.5 / 

7 = (t o /p) 7/v 

Blasius' empirical shear stress relationship is appropriate. 


T = 0.0255 pu 2 (v/u*f ) 0,25 


(3-4) 


where - R for pipe flew 


- 6 for boundary layers 


The boundary layer thickness implied by (3-4) and a 1/ 7 power-low profile is 


6 = 0.3 76 X/R 


0.2 


(3-5) 


From equation (3-3) 


*(^) = T q /p = (u*) 2 /v 


(3-6) 


In terms of real distance from the wall, equation (3-1) represents most of the boundary 
layer, therefore the following computation procedure is suggested. A fictitous wall is 
assumed to be 0.0005 feet away from the real wall, and no flow is assumed to occur 
between the two walls. Equation (3-1) is valid at .0005 feet from the wallf hence, if u is 
calculated with a slip boundary condition, u* is determined. Equation (3-6) is used to 
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calculate the velocity gradient at the wall. Since equation (3-1) is cot explicit In u*, the 
approximation 


» 

u 


0.1662529 u 


0.867325 


/(y/v) 


0. 132675 


(3-7) 


ia used. These equations determine the velocity gradient and shear stress at the walL 
An eddy viscosity is UBed to determine both the local shear stress and the variation of 
this stress with distance from the wall, y. 


U T = 0.07 u *£. (FR) + y (3-8) 

where 

FR = (y/0.3S' for 0< (y/fi, ) <0.3 

FR = 1. 0 for (y/£ ) >0.3 

Equations (3-7) and (3-8) and the momentum equations were used to calculate the 
turbulent boundary layer over a flat plate between 1 and 2 ft running length over the 
plate. A turbulent boundary layer was assumed at the leading edge of the plate. The 
flow was air with a free steam velocity of 100 fps (this is an approximate Reynolds 
Number of 10^). By adjusting the constants in equations (3-7) and (3-8), the profile at the 
end of the plate was predicted to be that shown in Fig. 3-2. The fit is very good, 
especially near the wall; the calculated wall shear stress is within 5% of the correct 
answer. This procedure 1 b accurate enough to extend its development to more 
geometrically complex flows. The reasons for the necessity of adjusting the constants in 
equations (3-7) and (3-8) and for the lack of better fit at ^s near the free stream side of 
the boundary layer are still under investigation. 

Research to continue these analyses until cylinders in cross-flow can be accurately 
simulated is in progress. 
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Fig. 3-2 Turbulent Boundary Layer Calculation 
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4. HGM FLOW ANALYSIS 


4.1 Introduction 

The analysis of the flow environment surrounding the SSME LOX posts requires a 
definition of the flow field in the HPFTP transfer tubes exit planes. The exit plane flow 
is development by defining the flow conditions immediately downstream of the turbine 
and computing the flow field through the turnaround duct, hot gas manifold and transfer 
tube. This section discusses the computation of the flow field in the turnaround duct, hot 
gas manifold and transfer duct for a two-duct configuration. 


4.2 Configuration 

The configuration analyzed consisted of the FMOF turnaround duct, the "Phase 3" two- 
duct hot gas manifold and the "Version B" transfer tube which includes the flow 
separator. The effects of turbine-induced swirl were neglected at the direction of the 
customer? hence, a plane of symmetry between the 2 transfer tubes was incorporated. 


4.3 Flow Conditions at Inlet 

The flow conditions at the inlet to the turnaround duct were specified by the customer. 
The fluid in the system was air at 530°R flowing at 72 lbm/sec. The pressure across the 
inlet was described by the equation. 


f 2 ^ 

P = 190.0 0.98 + 0.0441 Sin (-) 

l 2 J 


ps l a 


(4-1) 


where $ is the angular location which ranges from 0° between the transfer tubes and 
180° on the plane of symmetry on the side fartherest from the transfer tubes. The 
velocity profile in the TAD inlet is defined by the equation 
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V * V{ 1.0 + 0,04 Co#4>) {4-2) 

where V is the average velocity of any angle $ and V is the average velocity over the 
entire xnlct. The velocity has no cross flow component due to the assumption of no 
turbine-induced swirL The turbulent viscosity was specified as 10,000 times the 
molecular viscosity for air. 

4.4 Computational Grid 

The configuration described in subsection 4.2 was modeled using 10, 724 nodal grid 
points. The grid points depicting ths boundaries of the configuration are shown in Figs. 
4-1, 4-2, The struts and posts in the turnaround duct have been darkened in to clarify 
their locations. The computational grid in the plane of symmetry at the 0° (between the 
transfer tubes) and 180° (far side) positions are shown in Fig. 4-3. The inlet to the 
turnaround duct has been artifically moved upstream to avoid influencing the flow in the 
180° bend by the prescribed inlet flow conditions. Figures 4-1 through 4-3 illustrate that 
all of the salient features of the configuration have been incorporated into the grid. 
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Fig. 4-2 View of Grid From Inside 



COMPUTATIONAL GRID AT 160 DEG. PLANE OF SYMMETRY 



ACTUAL INLET 


COMPUTATIONAL GRID AT 0 DEG. PLANE OF SYMMETRY 


ACTUAL INLET 

/ 



FLOW SEPARATOR 


Fig. 4-3 Grids in Plane of Symmetry 
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4.5 Initial and Boundary Condition* 

The inlet conditions prescribed in subsection 4.3 were applied to the artificially displaced 
inlet shown in Fig. 4-3, The viscosity of 10,000 times the molecular viscosity of air 
results in a Reynolds Number of 300 to 400 and, hence, laminar flow. Therefore a 
laminar parabolic velocity profile was superimposed on the average velocity distributed 
defined by Equation . i. 

No-slip boundary conditions were specified for all solid wallB and tongency, or free-Blip, 
boundary conditions were applied to the plane of symmetry. The mass flow rate in the 
exit plane of the transfer tube was held fixed at 72 Ibm/sec. The total conditions at the 
inlet were held fixed, thereby allowing spurious signals to pass upstream and out of the 
problem. 


4 , 6 Results 

The flowfield for the TAD, HGM and transfer tube described above was computed using 
Continuum’s VAST code. The problem required 11,000 time steps before a converged 
solution was obtained. The results of the study were presented in detail to the customer 
on November 28, 1984. A summary of the results will be discussed in this section. 

The total pressure drop through the system is presented in Fig. 4-4 and shows a drop of 18 
psi in the 180 bend of the turnaround duct, a 14 psi drop through the struts, and a total 
drop through the system of 48 psi. Static pressure drop through the system was about 35 
psi. The pressure variation in the circumferencial direction in the HGM bowl inlet is 
shown in Fig. 4-5 and indicates a variation of 30 psi. Exactly how much of this result is 
affected by the assumed inlet pressure variation (8.38 pBi) is unknown but it appears to be 
small. 

The pressure distribution in the cross section in the plane of symmetry between the 
transfer tubes is presented in contour form in Fig. 4-6. The figure shows significant 
pressure drops in the 180° bend of the TAD and through the strut?.. The pressure 
variation in the bowl is small except in localized areas. 
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The pressure in the exit plane of the transfer tube varies only about 4 psl, hence, the exit 
plane pressure distribution is not presented in this discussion. Instead, velocity contours 
in the transfer tube are presented in Fig. 4-7 to show the nonuniform velocity 
distribution. A small area of high speed flow appears in the outer and upper portion of 
the tube. 

Figures 4-4 through 4-7 show large pressure gradients in the turnaround duct and email 
pressure gradients in the transfer tube indicating an improved design over the 3 duct 
system currently in use. 


5. CLOSURE 

The analyses of single and multiple LOX posts will be continued. The laminar flow 
analysis of the TAD/HGM/transfer tubes was completed; the turbulent flow case will be 
analyzed in the next reporting period. 
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ocity Contours Midway In Transfer Tube 
(Looking into Tube from LOX Posts) 


